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Abstract 

This work is intended as an attempt to study the non-perturbative renormahzation of 
bound state problem of finitely many Dirac-delta interactions on Riemannian manifolds, 
S'^ , , and . We formulate the problem in terms of a finite dimensional matrix, called 
the characteristic matrix The bound state energies can be found from the characteristic 
equation <!>(— 1/^)^4 = 0. The characteristic matrix can be found after a regularization and 
renormahzation by using a sharp cut-off in the eigenvalue spectrum of the Laplacian, 
as it is done in the flat space, or using the heat kernel method. These two approaches 
are equivalent in the case of compact manifolds. The heat kernel method has a general 
advantage to find lower bounds on the spectrum even for compact manifolds as shown 
in the case of §2. The heat kernels for if, and are known explicitly, thus we can 
calculate the characteristic matrix Using the result, we give lower bound estimates of 
the discrete spectrum. 



1 Introduction 

It is well-known that the exactly solvable Dirac-delta interactions on the plane and 3- dimensional 
Euclidean space in quantum mechanics give rise to some unphysical results for physical observ- 
ables, i.e., bound state energy and scattering cross section are infinite and the problem is said 
to be ultraviolet divergent. Nevertheless, there is a systematic way to dispense with these in- 
finities by means of a so-called regularization and renormahzation, which is first introduced in 
quantum field theory for the same reason. This problem constitutes an analytical example of 
regularization and renormahzation in quantum mechanics so that it helps us to understand and 
deal with it in a more elementary context rather than field theory and it has been studied in 
the literature from several point of views PP - Moreover, a single point interaction in two 
dimensional flat space is an instructive example of dimensional transmutation in non-relativistic 
quantum mechanics El El El- That is, the original Hamiltonian does not contain any in- 
trinsic energy scale due to the dimensionless coupling constant in natural units. Nevertheless, 
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a new parameter fj,"^, which species the bound state energy, must be introduced after the renor- 
mahzation procedure which then fixes the energy scale of the system. (A detailed discussion of 
dimensional transmutation in nonrelativistic quantum mechanics is given in a relatively recent 
article [IE])- 

In this study, we consider a bound state problem in which a non-relativistic particle living in 
a Riemannian manifold (in particular H^, and H^) interacts with finitely many Dirac-delta 
interactions. Similar to the corresponding bound state problem on and M^, we encounter 
divergences in this case as well. The main purpose of this paper is to show how to non- 
perturbatively regularize and renormalize the problem by means of heat kernel (even in the 
case where we do not have an explicit expression for it). After the renormalization, we estimate 
a lower bound for the ground state energy for each particular Riemannian manifold. This 
problem on two dimensional Riemannian manifolds, such as §^ and also displays a kind 
of dimensional transmutation ^Hj, where new energy scales different from the intrinsic energy 
scales of the system appear after the renormalization. We will briefly discuss it in sections 13.11 
audio 

Many body version of this problem on and is known as the formal non-relativistic limit 
of the scalar fleld theory in (2+1) and (3+1) dimensions. All these are extensively discussed 
in ^7]. Our primary motivation here is coming from the question how the renormalization 
method for the singular interactions in quantum mechanics would be performed on Riemannian 
manifolds, hoping that this may help us to understand the problem in the realm of quantum 
field theory. However, we shall postpone the discussion of the many body extension of it for 
future work and study first the one-particle Schrodinger problem. 

The paper is organized as follows. In section |2l we first define the bound state problem on 
compact and connected Riemannian manifolds and reformulate the problem in terms of a finite 
dimensional matrix $, which we will call as the characteristic matrix ^7j. Then, we emphasize 
the relation of the characteristic matrix with the corresponding spectral functions, resolvent 
and heat kernel. This allows us to reformulate the renormalization in terms of heat kernel. 
After that we continue to the discussion in the following sections by working out concrete 
examples. In section El we consider the delta interaction problem on §^ as an example for 
compact and connected manifolds. Considering the properties of the operator $ and using 
some properties and upper bound estimates of the heat kernel, Gersgorin theorem allows us to 
estimate a lower bound for the ground state energy of the system. In section |3J we apply the 
similar methodology, developed in the section of heat kernel method for S^, to the non-compact 
manifolds, such as and and show that the methods developed for compact manifolds work 
for some particular non-compact manifolds as well. Therefore, we renormalize the problem on 
hyperbolic spaces and give estimates on the ground state energy of each system. 

2 Renormalization of Finitely Many Dirac-Delta Inter- 
actions on Compact and Connected Riemannian Man- 
ifolds (M, g) 

The canonical quantization on non-trivial manifolds is known to have some ambiguities in 
quantum mechanics. For the path integral approach to the quantum system, the ambiguity in 
the canonical formalism is replaced by the undetermined parameter A and it can take various 
possible values ^H]- We remove this term for simplicity in all our examples, in which the 
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curvature term is constant and it corresponds to an overall shift in energy levels so that we can 
safely set A to be zero. 

Now, we consider a non-relativistic point particle living on a Riemannian manifold M in- 
teracting with a finite number of delta interactions located on the manifold and study bound 
states of the problem. We first investigate the delta interactions on a compact and connected 
Riemannian manifold (M, g) without boundary, of dimension D = 2,2) with the Riemannian 
metric g. The kinetic energy operator on Riemannian manifold (M, g) is just the Laplace- 
Beltrami operator or simply Laplacian, which is defined, in local coordinates 
for a neighborhood in the manifold, as follows: 

where g^fi is the metric tensor and g = [gap)- We shall usually denote the Laplacian as Ag to 
specify which metric structure on Riemannian manifold it is associated with. 

The spectral theorem 120] states that the eigenvalue problem Ag(f)i = Xi(j)i on a com- 
pact and connected Riemannian manifold (M , g) has a complete orthonormal system of C°° 
eigenfunctions 00, 01, • • • in L'^{M) and the spectrum Spec (A^) = Spec (M, g) = {A/} = {0 = 
Ao < Ai < A2 < . . . }, with A; tending to infinity as / — >■ 00. As a corollary of this theorem, the 
Laplacian on (M, g) provides us with all the tools of Fourier analysis, so that we can expand any 
"sufficiently good" function ip{x) on M in terms of the complete orthonormal eigenfunctions 

i;{x)=^QM^), (2) 



with the normalization 



(j)i{x)(j)i,{x) >/det g dx^ A ... A dx^ = 5iv 



M 



where C^'s are expansion coefficients. Note that extra labels in the eigenfunction expansion 
must be taken into account if the problem admits degeneracy. Delta functions on M can also 
assumed to be represented by these eigenfunctions 

5^{x - Oi) = ^ (3) 

1>Q 

with tti & M and 6^{x — a^) being the D - dimensional normalized delta function at point Oj, 

/ 6^{x - tti) ^/detg dx^ A ... A dx^ = 1 . 
Jm 

A typical Hamiltonian operator in quantum theory consists of a kinetic term, the Laplacian A^ 
with the factor fv^ /2m, and a potential function of position, attractive delta interactions in our 
problem. The time-independent Schrodinger equation on M for the bound states of a particle 
under the influence of N attractive delta interactions reads 



,2 N 



2m 

i=l 



ip{x) = -v'^ip{x), (4) 



where gi e M"*" is the strength of the delta interaction at and —u"^ is the bound state energy 
of the system. If we substitute (0) and © into the Schrodinger equation, it yields 



l>0 



n 



N 



2m ^-^ 



i=l 



0, 



where Ai = tp{ai) for simplicity of notation. The fact that (pi's form a complete orthonormal 
system allows us to write Ci in terms of them: 



1 



N 



gi (f)i{ai). 



(5) 



1=1 



Substituting (jS]) into the definition of Ai 



l>0 2m 



Xi + iy' 



and grouping the Ai terms we find 
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\, _|_ 7/2 



A,, 



N 

E 

i=i L 



5(j ^ (f)i{ai)(f)l{aj 

9: ^ 



^ \, _|_ 7/2 



A, = . 



The observation that the preceding equation is linear in Ai permits us to write it naturally as 
a matrix equation 



where $( 



$(-z/^)A = , 
is called the characteristic matrix and defined as 



(6) 



2JV 



-V 



A, + 



/>0 2m ' 



0i(aj)0;*(a 



if 2 



if z ^ J . 



(7) 



As we shall see below that the resolvent is intimately related to it and this allows us to state 
that the equation det$(— i/^) = gives the bound state energies of our problem. In other 
words, this equation is considered to be the determining equation of the ground state energy. 
Unfortunately, this nontrivial eigenvalue problem can not be solved analytically, that is, we can 
not obtain an exact expression for the bound state energy for arbitrary since the characteristic 
matrix depends nonlinearly on the bound state energy. Indeed, the problem is even worse than 
that, because we have not a finite expressions in the matrix elements of ^ij{—u'^). Fortunately, 
there exist a way to redefine the problem so that the physical observables yield finite values 
with the help of regularization and renormalization. Before introducing this procedure for our 
problem, it would be good to review first the problem in flat spaces. The infinite sums in the 



characteristic matrix on ]R2 or M3 is then replaced by integrals. The idea in that case is to take 



Fourier transform of the wave function 



^lj{k)e 



ik.x 



(2vr) 



D 



and substitute into the Schrodinger equation. Then we find that the diagonal part of the 
characteristic matrix is 

1 1 f d^k 



where D = 2, 3. This integral does not converge as it stands. The well-known method to 
remove the divergence is to put a cut-off A to the integral's upper limit and consider the 
equation as a determining equation of bound state energy for a given coupling constant g. If 
this regularization is performed, we realize that as the cut-off goes to infinity, ground state 
energy becomes divergent. In order to get a physically acceptable result, one assumes that the 
coupling constant depends on this cut-off and performs the limit A — > oo in such a way that 
bound state energy remains finite. These infinities should be removed properly since all the 
physical observables are measured experimentally as finite quantities. The cut-off dependence 
of the coupling constant is chosen as 



d^k 

17'. 

|fc|<A 



^,(A) (27r)^ L|<^P + ;xf ■ 



The determination of this coupling constant is called renormalization. Now, we follow the same 
idea to remove the divergence from our problem. By using Weyl's asymptotic formula |21j . 
one expects that the diagonal term ^^>o If'^"'"*^ ^ in the above matrix does not converge and 

this will be explicitly seen for a particular manifold S^. For a general compact manifold, we 
introduce cut-off to the upper bound of the infinite sum and choose the coupling constant as 

where —fif is the measured binding energy to a single delta interaction. Then, we take the 
limit A ^ oo 



lim 

A^oo 



(10) 



and this should give us a finite result in two and three dimensions. Hence, the divergence has 
been removed and bound state energy becomes finite. A rigorous proof of this is not trivial, so 
we will stay at a heuristic level and study special cases only. 

As we will show in the next subsection, the heat kernel is intimately related to the charac- 
teristic matrix $ and this relation helps us to see easily which part of the matrix is divergent 
or convergent and then how to renormalize the problem non-perturbatively. Furthermore, heat 
kernel is especially very helpful to remove the divergences for our problem on non-compact 
manifolds, as we shall discuss in section |3J We will see that the above method can easily be 
extended to find the renormalized resolvent of the singular Hamiltonian. 



2.1 The Relation of Matrix $ with Heat Kernel and Resolvent 

The resolvent (or Green's function) and heat kernel play very essential role in establishing the 
connection between spectral properties of the operator and corresponding geometrical notions. 
Up to now, we have been dealing with a matrix $, and do not refer to resolvent and heat kernel. 
In order to see the relation between the matrix $ and heat kernel we consider the separable 



Hamiltonians H = Hq — Xl/li fi'il/j) (/«!' "where is a particular Dirac ket. We work out the 
resolvent formula of H in terms of Hq and assume that the two Dirac kets and \x) are 
related in such a way that the equality {H — z)\ip) = \x) is satisfied. Then, we have 



N 



\^) = \x) 



assuming complex number z ^ Spec(i/o)- Acting the operator {Hq — z) ^ on both sides and 
projecting it onto we obtain 



N 



Y^^,,{z)U,\^) = gl'{U\{HQ-zr'\x) 



where we define a matrix ^iAz) as ^ 



g7^-{f^\{HQ-zr'\U) lil=3 

-ii{n\iHQ-zr'\f,) iu^j. 

9i 



(12) 



After a little algebra, it is evident that 



(H - zy' = {Hq - zy' + (Hq - zy' 



N 



.i,j=^ 



(Hq - zy' 



(13) 



as long as ^ij{z)~^ exists. Such formulae were extensively discussed in problems associated 
with self-adjoint extensions of operators, notably by Krein and his school, and also for such 
singular interactions in flat spaces Therefore, our problem can also be considered as 

a kind of self-adjoint extension of the free Hamiltonian. It is defined through regulating (or 
controlling) the behavior of the wave function in the vicinity of these interaction points. 

If we take the matrix element of (jl3|) by projecting on to the Dirac kets {x\ and \y), we have 
found the resolvent kernel R{x,y\z) = {x\{H — zy^\y) corresponding to (fTT|) 



RixM^) = RoixM^) + JdxWRoMz: 

= RQ{x,y\z) + ^ / dx'RQ{x,x'\z)fi{x') 



N 



5^ /,(x')$.,(.)-V,(z/' 

.«J=1 



Ro{y\y\z) 
dy'Ro{y',y\z)fj{y') 



By choosing the functions /i(a;)'s as bump functions centered at x = at such that the sequences 
of the functions admit the limit fi{x) — > S^{x — ai) (in the appropriate topology), it turns out 
that 
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R{x,y\z) = RQ{x,y\z) + ^ RQ{x,ai\z) ^ijiz) ^i?o(aj,2/| 



z . 



(14) 



^ There is no confusion in notation because we will see that this matrix $ is exactly the same matrix considered 
in the previous sections. 



The important point to note here is the relation between the resolvent operator, defined on an 
infinite dimensional space and the characteristic matrix, defined on a finite dimensional space. 
This allows us to find the bound state spectrum of the separable Hamiltonian operator H with 
the help of a finite dimensional matrix <^{z)~^. Since discrete spectrum is the set of complex 
numbers such that the resolvent does not exist, this proves that the equation det $ = gives 
the bound state spectrum of our system. The fact that the free Hamiltonian is bounded from 
below allows us to write the free resolvent operator as an integral for 3^(2) < 



{Ho-z)-' = ^l e"^(^^^"")cit, (15) 



the result of which should be continued analytically to its largest set in the entire complex 
plane. As a consequence of this, the free resolvent kernel is 

Taking the limit fi{x) 6^ {x — ai) , it results in 

1 zt 

Ro{ai,aj\z) = {ai \ {Hq - \aj) = t Kt{ai, aj) dt , (16) 

" Jo 



where Kt{ai, aj) is the so-called heat kernel, and the operator e ^'»Jv2™ "y is the formal solution 
to the heat equation ^^llin]- Hence, the matrix $ is written in terms of the heat kernel in the 
following way, 

{1 t 
- - / e^Ktitti, Qi) dt Hi = j 
o, 1 A ^ (17) 
— - - er^ Kt{ai,aj)dt if z ^ j. 
9i JO 

The matrix is exactly the same matrix mentioned in the previous sections. This can be 
shown easily from the spectral theorem for compact manifolds 

Kt{a,,a,) = ^e-^"'[^]0z(aO0;(a,-), (18) 
which converges uniformly on M x M for each t > 0: 

rco 7 . pco / \i t] rit 



z>o 



(j)i{ai)(j)^{aj 



where summation and integral are interchanged since summation converges uniformly. This is 
the same result for z = —v"^ that we already obtained for non-diagonal part of the characteristic 
matrix in the section |21 One can understand how the non-diagonal part of it in ((7j) is convergent 
by using the smooth behaviour of the heat kernel and the integral e"^ Kt{ai,aj) is convergent 
for Qi 7^ aj. However, the asymptotic behaviour of the heat kernel as t ^ 0"*" for every point x 
on a compact manifold M jT^] is given by 

°^ fht\' 

(19) 



^ ^ k=0 ^ 
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where D is the dimension of the manifold and the Uk{x^x) are functions given in terms of 
the curvature tensor of M and its covariant derivatives at the point x. This result shows 
that diagonal part of the heat kernel as t — *• 0+ for D = 2, 3 leads to a divergence since 
uq{x,x) = 1 (there is no infinities for D = 1 as it can be easily realized). In other words, 
the sum in the diagonal term in $ is divergent while the sum in the non-diagonal term is 
convergent. However, we have already shown that bound state energies are related to the 
characteristic matrix, i.e., det^{z) = contains information about bound states. If some of 
the elements of the characteristic matrix have infinities, it is impossible to get sensible bound 
state energies for our problem. Before establishing the renormalization of our problem with the 
help of heat kernel, we must indicate why this problem occurs. Although the delta interactions 
may approximately describe a system in which a particle interacting with a point-like centers 
when its de Broglie wavelength is large compared to the typical range of a potential, we have 
not encountered in nature this type of contact interaction. This means that the substituting 
the Dirac-Delta interactions into the Hamiltonian for D = 2,3 directly is not a proper way. 
Therefore, we must modify our problem such that it has a finite range and then consider the 
zero range limit. In our renormalization method with heat kernel, short range is replaced with 
the short time as we will see. 

We introduce a small constant e, in the lower limit of the integral. We then take the limit as 
the cut-off e goes to zero in such a way that the experimentally measured ground state energy 
remains finite. This requires that some quantities in the problem, e.g. coupling constant, should 
have a cut-off dependence in a definite way. For our problem, we naturally choose 

Qi {e,fii) = - / e ft Kt{ai,ai)dt . (20) 



After performing the limit e ^ 0, we have the renormalized characteristic matrix 



t 

oo 



e ft —eft 



dt if i = i 

(21) 



e n Kt{ai, ttj) dt lii ^ j 







where '^{z) < and ^ij{z) can be analytically continued to its largest set in the entire com- 
plex plane. One can naturally ask whether the renormalization performed with heat kernel is 
compatible with the one introduced in section |21 The answer is affirmative and one can easily 
show that the cut-off A for the infinite sum introduced in section |21 corresponds to the cut-off e 
for the lower bound of integral in the heat kernel method. This can be realized easily by using 
the spectral theorem in the diagonal part of equation (PT|) and taking z = —v^: 



z>o 

=^7-l(e,/i,)--5^0K«.)0r(«^) / e^e-^^A^dt 
^ i>o 



(22) 



where we have used the uniform convergence of the sum. Now, in order to remove the divergence, 
we can naturally choose the coupling constant as 

-| poo _ 2^ 2 r n 

9i\e,^,) = -Y^<Pi{ai)<Pl{ai) e^e-^'^ii\dt. (23) 

7\n ^ € 
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Then, we have 



hm , ^ 



^(f)i{ai)(f)^{ai] 



-j4t_ 



dt 



--y 



(24) 



which is the same result we would have obtained by the eigenfunction expansion by introducing 
a cut-off A (equation (|iup). After finding the renormalized characteristic matrix, the resolvent 
can be written explicitly 



N 



R{x,y\z) = Ro{x,y\z) + ^ Ro{x,ai\z) ^ij{z) ^Ro{aj,y\z) 



where 



Ro{x,y\z) 



h 



e'r^Kt{x,y) dt . 



(25) 



(26) 



Once we have given the resolvent of an operator, all the information about the operator is 
contained in it. Nevertheless, it is instructive to check that the wave functions can also be 
obtained and they are normalizable. We write the normalized wave function with a cut-off A 
and then take the limit A — > cxo, this way we will not get a vanishing wave function. So the 
normalization constant can be found easily 



N „ A 

i,j=l l,l'=0 



(f)i{ai)^^{x) ^^,{aj)(l)iix) 



N 



5:<,.(A)a(A)4-(A)A(A)E#r^ 



(27) 



One expects from the Weyl asymptotic formula that the wave function is not normalizable if 
we are on a space of dimension bigger than three. Moreover, we can see that the summation 
over the eigenmodes is exactly the derivative of $(— z/^) with respect to z/, hence we get: 



ic(A)r 



(2^ 



Performing the limit A ^ cxd, the properly normalized wave function of n*^ state becomes 



N 



9$. 



r,s=l 



N 



0;(a,)0;(x) 



i>o i=i ih + ^l) 

where z/„ is the n*'* root of the energy equation det $(— z/^) = 0. This can further be simplified 



to an expression in terms of the heat kernel 



N 



r,s=l 



du 



N 



X 



N 



r,s=l 



l>0 i=l 






_ 1 




2 />oo 




/ 




Jo 



h 



poo 

■^0 i=l 



dt 



(29) 



in which one can easily see that tpnix) is finite. 

3 Finitely Many Dirac-Delta Interactions on 

Since the simplest and one of the most familiar compact manifolds is the sphere S^, we shall 
work out the problem of point interactions on a sphere as a concrete example. Suppose that 
point interactions are located at the points given by the local coordinates {9i, 4>i)f=i on a sphere 
of radius R. Then, the Schrodinger equation for the bound states of a particle living on the 
sphere under the infiuence of attractive delta interactions becomes 



2m 



N 



1=1 



ip = —V ip, 



(30) 



where A§2 is Laplacian on the sphere in spherical coordinates 

1 92 



A, 



1 d f . B 
smu- 



i?2 sin e 36 



09 / i?2 sin^ 6 dcP 



(31) 



e,) t5(^ 0,) ^j^^ dimensional delta function on the sphere centered 

R-^ sin 9 ^ 



and 5'^{9-ei,(j)- 

at {6i,(j)i). It is well known that spherical harmonics Y^^ are eigenfunctions of the Laplacian 
Ag2 with the eigenvalues /(/ + 1)/R^ and form a complete orthonormal basis on S^, order 
to be consistent with the standard normalization of spherical harmonics, we choose (pim 
From the following identity 

I 



R 



2/ + 1 

2/ + 1 
47r 



Pi (cos 9i cos 9j + cos(0,j — 0j) sin 9i sin 9j) 
d? 



Pi I 



2 



in 



where dij = -j^ 



fjl e [0, 2] being rescaled distance between point centers with radius of 



the sphere R, the matrix $j 



-u"^) in ((7j) becomes 
1 21 



47ri?2 



/(/ + l) + Z/2 



9j 



9r 47ri?2 



-R2 



i>0 2m/?2 
2/ 



1 



+ + 



1 



d^. 



(32) 



^ 7^ J- 



i>0 2m_R2' 

It follows easily from the Cauchy-MacLaurin integral test that the infinite sum 

21 + 1 



-y 



ft2 



l>0 2m_R2 



/(/ + 1) + Z/2 



is divergent. To get a sensible results for our problem, we must modify our original problem as 
outlined in section |^ Therefore, considering our problem in the light of this method, we first 
define the coupling constant gi as a function of the parameter A (cut-off). Then, by choosing 
g~^{Ays naturally 

1 A 2/ + 1 

where /ij is experimentally measured value of bound state energy for the single delta interaction 
and taking the limit A — > cxd of the difference, we have obtained 



lim 



1 ^ 



21 + 1 



1 ^ 



21 



4-^' + 1) 

V 



where = ^^^^ . The function here is defined as 



X G 



where if's are the harmonic numbers, commonly defined on integers as Hn = Ylt=i k '"^^ 
be extended by analytical continuation to its largest domain in the entire complex plane as 



iIj{z + 1) + 7, where iIj^z) 



T{z) 



being the digamma function and 7 being the Euler- 



Mascheroni constant. The digamma function has several useful integral representations 
some of which are 



-zt 



\ogz + 



t 1 - e-* 
f 00 / 



dt 







1 

- - 1 
t 



-zt 



dt 



(33) 
(34) 



where 'R{z) > and these can be useful for the estimates of its upper and lower bounds. Due 
to the Schwarz reflection principle of harmonic numbers {Hz = H^), the function (f){x) is real 
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valued 



for all X G 



It is also easy to check lim 



1 in the non-diagonal 



3. (A) 
A^°o 3, (A) 

part of ()H2|1 , simply because of their same form of the divergence. Then, the renormalized 
matrix ^f— 



for bound states can be eventually written as 



V 



i>0 



2/ + 1 



(35) 



^ 7^ J 



By the analytical continuation of the characteristic matrix to its largest domain in the entire 
complex plane, we have 



21 



P, 1 



d^. 



J 



(36) 



from which we can write the resolvent equation (jl4p . Hence, we have obtained a well-defined 
formulation of our problem, that is, the infinities have been removed. Moreover, we see that the 
problem realizes a generalized dimensional transmutation. In this case, the coupling constants 
gi have the same dimension as ^ by dimensional analysis. In contrast to the fiat case, we 
have one more parameter R coming from the geometry of the space. Thus, we expect that 
the system must have an intrinsic energy scale ^^^2 as well as terms. However, after the 

renormalization, we obtain a set of new dimensional parameters /if. Hence, the first set of 
scales we expect by naive dimensional analysis at the beginning is not sufficient. Instead, a 
specific combination of all these parameters together determine the scale of our problem. This 
means that delta potentials on a sphere is an example of a kind of dimensional transmutation. 
However, there is a slight difference, especially in the case of single delta attractor: in the 
fiat case there is no combination of dimensional parameters to come up with an energy scale, 
whereas in the case of a sphere we have a geometric length scale R which already defines an 
energy scale • The dimensional transmutation is most striking in such cases where there 
is no intrinsic energy scale. 

In order to estimate the non-diagonal part of the matrix $ for sphere S^, we follow a different 
strategy, using the heat kernel. 



3.1 Lower Bound of Egr by Heat Kernel Method for 

Heat kernel Kt(x,y) is the unique fundamental solution to the heat equation -^Agcf) = —fvpt- 
It has the symmetry {Kt{x,y) = Kt{y,x)) and semi-group property [121 12n]- As well as being 
a useful computational tool in establishing the existence and some of the properties of the 
spectrum of the Laplacian of the eigenfunctions on a Riemannian manifolds, it is very helpful 
to understand the nature of the divergences for our purposes, as we have shown in the previous 
section. 

By means of the relation !^'21\i and explicit form of the heat kernel, one can calculate the 
matrix However, there are some situations in which one can not calculate the heat kernel 
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explicitly, e.g. we do not have an explicit expression of the heat kernel for two dimensional 
sphere. In this case, one can still find some bound estimates on matrix $jj without having 
explicit form of the heat kernel, instead some properties of it. In order to analyze this for S^, 
we will use some estimates on the heat kernel, based on a work by Li and Yau [21]. Let us 
recall the corollary of the theorem (3.1) in 

Let M be a complete manifold without boundary. If the Ricci curvature of M is bounded 
from below by —K, for some constant K > 0, then for 1 < a < 2 and < e < 1, the heat 
kernel satisfies 



Kt{x,y) < 



C{e) 



V{x,Vt)V{y,Vt) 



where V{x,r) = fi{B{x,r)), B{x,r) is the geodesic ball of radius r centered at x G M and 
d{x, y) is the geodesic distance between two points x and y on the manifold. The constant C7 
depends only on the dimension of the manifold D, while C{e) depends on e with C{e) ^00 as 
e ^ 0. When K = 0, the above estimate, after letting a — > 1, can be written as 



Kt{x,y) < 



C{s) 



V{x, ^rt)V{y, Vi) 



(37) 



Since §^ satisfies the conditions above as a particular case, this corollary can be applied to 
it as well. On the other hand, we have a different purpose from the original corollary of the 
theorem for the estimates on the upper bound of the heat kernel, in which the sharp estimate 
for the heat kernel is found. Instead, we are trying to find a best lower bound of the ground 
state energy of the system. Therefore, we shall modify the original corollary in (2^. Using this 
theorem with relaxed condition < e < 1, we have found the upper bound estimate for heat 
kernel of sphere in our problem: 



< 



C'i6) 



»L 

g D{S)ht 



V X 



(3^ 



where 



and 



C\6) = {l + 6)^exp 



1 



+ 



1 



1 

45(1 + (5) (1 + 2(5) ' 25(2 + 5) ^ 45 



(39) 



D(5) = 4(l + 25)(l + 5)^ (40) 

5 is merely required to be positive. When we want to find a lower bound for the energy, the 
numerical values of the coefficients C"(5) and D{6) will be determined explicitly. It is easy 

to see that V (x, = V {v^ = ^^R' 



cos 



m 

2m/ 



^ j as long as < t < 



2r2 



h 



For t > 2S2I^J^^ we have V ^x, sj^^ = V (y, \f^^ = AttR'^. According to our corollary and 
positive definiteness of heat kernel, the following integral has an upper bound: 



"ft Kt{ai,a.j)dt 



< 



C'i6) 

h 



g D{S)ht 



2 2 
ft 



27ri?2 1 - cos 



ht 

2mW 



dt + 



C'jS) 



i] V t 

g D(S)ht ft 



dt 



where we have taken z as —v^. With the help of the identity 1 — cos y 2mB? 
and the inequahty ^t^^ < ^ for < ^ < 7r/2, we obtain 



ht 



2sin^ 



ht 



8mR2 



-i^r. / X , rmrC'iS) 
e » Kt{ai, aj) at < 



2^2 



g D(S)ht h 
t 



dt 



C'jS) 



e "^H^Ifit ft dt 



Evaluating these integrals, we find 

poo 



-V 



h 



e ft 



Kt{ai,aj)dt < C'{5) 



— Koia^.u) + - 



where 



a,: 



8m 



(41) 



and Ko{x), Ki{x) are modified Bessel functions. This shows us that the infinite series in the 
non-diagonal part of the characteristic matrix is finite and bounded from above according to 
(0^. In order to find a lower bound for the diagonal part, denoted by V, of the matrix $ for 
sphere S^, we first recall how the diagonal part of the matrix $ appears in : 



E 



2/ + 1 



1=0 2mB? 



E 



21 + 1 



1=0 2mB? 



/(/ + !) 



> . 



Instead of calculating explicitly this limit as we have done in section El we estimate a lower 
bound of it by means of integrals replaced by the sums as follows 



r lim 

47ri?2 A^oo 



A+1 



2t + 1 



n? 

2mW 



t{t 



dt 



After taking the limit we find 



v.: > 



m 



log(z///ii 



2t + 1 



?i2 



2m 



^t(t+l) + Z/2 



dt 



and using the estimate for logarithmic functions in [22] 



X — 1 

log X > for a; > , a; 7^ 1 



(42) 



X 



we obtain 



m 



log(z///ii 



47ri?V2 



> 



m 



TT 



AttR'^u^ 



> 



(43) 



(44) 



For positive definiteness, we have assumed u is sufficiently large, which is not a particularly 
restrictive condition. In fact, one can try to find sharper estimates by means of the integral 
representations of digamma functions ()33|) and ()34|) without this assumption. However, the 
estimated functions in this case are too complicated to suggest a bound for ground state energy. 
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A well-known theorem in matrix analysis, called Gersgorin Theorem j2n] states that all the 
eigenvalues Aj of the renormalized matrix $ are located in the union of discs 



N 



- <^u\ < Rim ^ Gi<^) , 



(45) 



1=1 



where R'ii^) = X]i^j=i \ 1 < i < N. If we want A = not to be an eigenvalue, then 

none of the discs should contain A = 0. Then, we should impose 



TV 



V,{u)\>Y,\JC^j{u)\ 



(46) 



for all i. This is possible for a critical value v > v* since the left hand side is an increasing 
function of v and the right hand side is a decreasing function of it for a given d and A^. In fact, 
this inequality obviously provide a lower bound for the bound state energy by just plotting the 
functions on both sides in spite of how complicated the form of functions are. However, we shall 
try to find an explicit expression for the lower bound of the ground state energy depending on 
the number of delta interactions. In order to achieve this, we choose v such that; 



> 



m 



mjii 



1 



N 



(47) 



where we have used the monotonic behavior of the functions in Vij and JCij and defined n = 
maxj fii and a = miuj^j aij or d = miuj^j dij. From the integral representations of the Bessel 
functions for z G M"*" 



Ko{z) 



-z cosh t 



dt 







oo 



z I e -'""'^^^sinh^t (it 
'o 

2t 



and using the inequalities y < cosht, sinh t < ^ for all t G M"*", we can find the upper bounds 
for the functions Kq and Ki 



Ki[au) < e 2 — 



(48) 



where ai> G M"*". Considering the estimated bounds for Bessel functions, it is easy to see that 



m 



{N-1)C'{6) 



e 2 /ii e 2 ^j^a 



+ 



> {N-1)C'{6) 



Using the argument v > v^/ij? in equation and last inequality, we impose the following 
inequality with the help of Gersgorin Theorem: 



u > n + 



2V2 



+ {N -l)C'{6)e-'^ 



a 



2V2 



(49) 



Let us make the following reasonable assumptions and take these for granted for the present 
(we will later show that they indeed satisfy these conditions by finding the extremum of ground 
state energy with respect to the parameter 6) 



(50) 
(51) 
(52) 





< 


27r2 


4 


a 


fJ'R 


< 


27r2 




a 


1 


> 


1 

5 ' 



so that the inequality becomes 

> ^ + 



2V2 



+ 37rVd(A^- l)C"(5)v/^D(5) e ^^^^ 



(53) 



from which we conclude that there exists a critical value u > u* for a given such that A 7^ 
and then, the ground state energy cannot be less than —u*'^: 



Egr >-u* 



fJ'R 

2^2 



+ 5/irfiy 



5 



(54) 



where W is the Lambert W-function, also called Omega function or product-log function 
Now, we choose 5 in such a way that the energy bound take its minimum value. This is 
accomplished if 5 is chosen approximately 0.508, which is independent of the parameters in 
the problem. This independence can be easily realized from the form of inequality (|53p . By 
substituting the values of C'{5) and -D(5), we estimate a lower bound for the ground state 
energy: 



Egr > —V 



*2 



287r^ (A^ - 1) e 



2-/2J 



(55) 



By using this value of 6 and the fact that d < 2tiR, the consistency of the assumption we made 
can be shown easily. Finally, we shall consider the large N behavior of the ground state energy. 
The asymptotic expansion of product-log function W |2H] for large z is given as 

(56) 



W{z) ~ log 2 — log log . 



Hence, this leads to 

Egr ~ -/i^ [log (iV) - log log {N)f . (57) 

The method we have introduced for the two dimensional sphere §^ can also be applied to 
a general compact manifold. The main idea is based on finding an upper and lower bound for 
the characteristic matrix or heat kernel (based on the work by Li and Yau). Then, Gersgorin 
theorem allows us to estimate a lower bound for the ground state energy. 



4 Finitely Many Dirac-Delta Interactions on Hyper- 
bolic Spaces 



The hyperbolic space M."" is defined as maximally symmetric and simply connected complete 
n-dimensional Riemannian manifold with a constant negative sectional curvature —1/R, which 
is also in some sense considered to be the negative curvature analog of the sphere E>"'. We shall 
deal with the delta interactions on the hyperbolic spaces and in the following sections. 
The method developed in the previous sections for §^ will be useful as well for the hyperbolic 
spaces. The heat kernel on hyperbolic spaces [22] ? written in terms of dimensionless quantities: 



(58) 



V2 



2 



^^^"'^^ = (4^^^"' L (cosh.- cosh rf)V2 ^ ' 

where d = dist(x, y), geodesic distance between two points x and y on H". 

Although spectral theorem and asymptotic expansion of heat kernel discussed in the previous 
sections may not be valid for general non-compact manifolds, we shall demonstrate that for 
the specific examples in non-compact manifolds, such as and H"^, our viewpoint still works. 
It would be desirable to show the equivalence between the eigenfunction expansion and the 
heat kernel method for the regularization in non-compact manifolds rigourously. Nevertheless, 
we have not been able to do this. The main idea is similar in spirit to the renormalization 
procedure introduced for the compact manifolds. 



4.1 Finitely Many Dirac-Delta Interactions on Hyperbolic Space 

In the hyperbohc space M.^ = {x & M.^\x3 > 0}, the geodesic distance d is defined as 

d{x,y) ^ , \x-y\^ 



cosh ■ 



R 



1 



2 X3 ys 



where R is the scaling parameter. The Schrodinger equation for the bound states of a particle 
living on If under the infiuence of N attractive delta interactions is 



2m 



N 



i=l 



(59) 



where Laplacian A^a in polar coordinates (x, 0, 0) 

1 (9^ 2coth^ d 



A 



H3 



R^d^P^ 



+ 



1 



R^ dip i?sinhV 



^§2 



(60) 



We have an explicit formula [SH] for the heat kernel of the three dimensional hyperbolic plane 
if written by using physical constants 



Ktix,y) 



R 



R 



exp 



ht 



2mR^ 



md{x, yY 
2ht 



(61) 
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such that as R 
kernel as 



oo, we can obtain the heat kernel on M'^. Hence we have the free resolvent 



(oil (Ho 



\aj) 



oo exp ( f 



M 



2mB? 



mS ■ 
a. 

2U 



dt 



"■1-3 

R 



exp 



1-4 



(62) 



where d 



d{ai,aj), fi% 



h? 2 



It follows easily that this term gives infinity 



when i = j, that is, the diagonal term in the characteristic matrix is divergent. Then, we can 
now proceed the regularization and renormalization schemes analogously for the hyperbolic 
spaces. However, the divergence in hyperbolic space is due to fact that the lower bound t 
of integral ()62j) is zero. Hence we regularize the divergent term by introducing a lower cut-off 
e, as we have shown in section ITTl we expect this should in some way related to the ultraviolet 
regularization. We next define the coupling constant as a function of this cut-off: 



^7 



lim 



1 



9i 



(47r)3/2;i|i?3 



du 



where the integration variable u = ^^^2 1 is introduced for simplicity. The natural choice for 
g-^{e) is simply 



9i 



1+^ 



du. 



where /Xj is an experimentally measured bound state energy for the single delta interaction and 
it helps us to keep track of the strength of point interactions. In e ^ 0"'' limit, we have found 
the explicit renormalized characteristic matrix for M.^ 



1 - 



1 + 



f^R sinh 



exp 



1-4 

f^R 



if i 



if i ^ j. 



(63) 



Then, we have the resolvent equation (fT^ with the free resolvent kernel Ro{x, y\z) for H'^ given 
by 

1 



Ro{x,y\z) = — 



d{x,y) 
R 



47r/i|i?3 sinh^^ 



exp 



I^R 



^^d{x,y) 



Z 

1- — 

^R 



(64) 



from which we can get all information about the system. Using the Gersgorin Theorem ()45j] 
for this matrix, and following the same ideas introduced for we obtain 



> {N 



d 



fiRR sinh ^ 



d 




where we have taken z = —u"^ and chosen d = miuj^j dij, and fi = maXjyUj. It turns out that 
this inequality indicates that there exist a critical u > u* for a given d and N for which this 
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inequality is satisfied and zero is not an eigenvalue. Therefore, the ground state energy cannot 
be less than —u*^: 



/i'-2/irf.//i2j + /i2 W 



"■^ V ^ d{N - 1) 
/2sinh| 



^'^ V ^ d{N 



i?sinh 



R 



(65) 



For the large N behavior of the ground state energy, the estimate becomes 

Egr 2/i,v/;4T7?[logiV - loglogiV] - fil [logiV - log log iV]' 



Now let us consider the two center case on the hyperbolic plane M.^ and assume again that 
their strengths (or bound state energies of each center) are the same. In this way, determining 
equation (det $ = 0) becomes 



'1 + 



jJ^R sinh I 



exp Wl + — 



If we expand it for small d we have 



f^R 



fJ'R 



1 - . 




Z/2 









from which we can conclude 



E. 



gr 



4^^^ 4 4 2 V /i| 



Similarly, for large values of d, the right hand side of the energy equation for two dirac delta 
interactions vanishes, so that we obtain the ground state energy E, 



gr 



4.2 Finitely Many Dirac-Delta Interactions on Hyperbolic Plane 

The geodesic distance on the hyperbolic plane is defined by 

cosh^ = 1+1^-^1' 



R 



2 ?/2 



where i? is a scale distance. Then, the Schrodinger equation for the bound states of a particle 
living on in the presence of attractive delta interactions is 



N 



— Ah2- V 5\e-e,,^- 

2m ^-^ 



■il) = —V ijj, 



(66) 



IP 



where the Laplacian A]fj2 in polar coordinates 0) is given by 

1 92 2coth^ d 1 92 



A 



H2 



The heat kernel for j^H with the proper physical parameters is 



(67) 



Kt{x,y) 



1— ft t „ 



re 4 1 t 



A / cosh r — cosh -'^^ 



(6^ 



One can check that this goes to the heat kernel oys.^ qjs, R ^ 00. Then, the free resolvent 
kernel is immediately obtained 



r. 1 



4:71 fi^R'^ ./fii rf: 



ft _t 



(ir. 



_r±2rnR^ 1 
re 4 ft t 



~R \l cosh r — cosh ^ 



dr 



dt 



w Y cosh r — cosh ^ 

We see that the diagonal term, which corresponds to dij = is divergent, as expected. Therefore 
we again repeat the similar regularization and renormalization procedure as we have done for 
H^. After introducing a cut-off to the lower limit of the integral 



^ii{z) = lim 

e— >0+ 



lim 



1 



9i lej 



9^ (ej 



47r/i|/22 Vcoshr - 1 
V2 



dr 



47r/i|j,i?2 sinhu 



du 



and by the natural choice for ^(e) 



9i (ej 



V2 



«. 1+ 



4Mf 



du , 

sinhw 



we have obtained the renormalized characteristic matrix for in the e — > limit 

1 



1 



V2 



'1 ifl 

' 2 



r, 1 



cosh r — cosh 



(ir 



if i = j 



if i 7^ J, 



(69) 



9n 



where ip is the digamma function. Then, we have the resolvent equation (jl4j) in which the free 
resolvent kernel RQ{x^y\z) for is given by 



Ro{x,y\z) 



r. 1 



dr. 



\ cosh r — cosh 



R 



(70) 



The integral on the right hand side is in fact one of the integral representation of the Legendre 
polynomials of second type P?] 



V2Qx ( cosh 



R 



-(A+|)r 



dr 



~R~ \ I cosh r — cosh 



d{x,y) 
R 



(71) 



which are defined for 3? (A) > —1 and in our case 3fJ (A) = 3? {^■\J^ ~ ~ ^ > — 1- Therefore, 
the free resolvent in terms of Q\ 



(72) 



Gersgorin theorem allows us to estimate the lower bound for the bound state energy as done 
for §^ and H'^. In order not to have zero as an eigenvalue, we must have 



V2 



2 



'1 



f^R 



2 



'1 



f^R 



dr 



coshr 



cosh ^ 



(73) 



for all i and we have taken z = —v^ and v > maxj /Xj. It is easy to see this inequality is satisfied 
for some values of u because the left hand side is an increasing function, whereas the right 
hand side is a decreasing function of u. However, it is not so easy to give an explicit estimate 
for i> by this inequality so we will estimate the functions on both sides. The inequality for the 
digamma functions [201 



ip{x) > logx X > 



X 



(74) 



which can be obtained from the integral representation ()34|1 . and ^ + yj + x'^>x for all x > 
helps us that we can find the following inequality by assuming u is sufficiently large 



2 



'1 



2 



4 /i| 



> 



log 



^^R 



I'r 









(75) 

Since the right hand side of equation (f7S|) is (A^ — l)y/2Qx ^cosh ^ j we can find an upper 

bound for this function, using another integral representation of the second type Legendre 
polynomials fT7\ : 



Qx { cosh ^ 



r(A + i) 



-t cosh 



R Kq it sinh 



di-i 



13 \ +\ 



R 



r dt 



(76) 
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where = and A = ly^l + ^ — |- Using the estimate for the function Kq given in 

equation (pSj). we obtain 



V2Q, 



cosh 



dij 



< 



2V2 



-*(™sh%- + ismh%-)^A-l 



V RJ r(A + l)sinh^^ io 
and the right hand side is just the Gamma function, then the estimate becomes 



V2Qx ( 



cosh ■ 



R 



< 



r(A + 1) (cosh ^ + I sinh ^] sinh 



R 



Using identity r(A + 1) = AFfA) and the assumption u/ fi^ > 1 and ^/l + ^ > — , we get 



V2Qx r 



cosh 



d 



< 



4^2 



V Rj (^-l)f 



^cosh ^ + I sinh ^ j ^ sinh ^ 

1^ 1 N ' 



< 



1/ 



sinh^ 



(77) 



Also, by choosing d = miuj^j dij and /i = maxj /ij, we easily find 



> 



log ^ o + '\/7 + ~ 



(7^ 



and 



logfcosh + | sinh -g^- 



e 



sinh ^ 



< 



sinh 



Therefore, we impose 



(79) 



1 



log ^ o + '\/7 + — 



> 



4(Ar - 1) e-^^-^)'°s(™^^^+^"°''^' 



sinh 



^0) 



It is immediately seen that there exists a critical value v > v* for a given c? and for which 
this inequality is satisfied and zero is not an eigenvalue. Last inequality can be written as 



— log 

\iR 



\ 



V 



\ 



> 



4(Ar- 1) 
sinh "I 



1^ 



-1 



If^>e V* '"rJ (independent of A^), then we have the lower bound of the ground state 
energy 

' 4(jV-l)Ae-'^/^ ''\ -, 2 
sinh - 



Egr > —y 



*2 



■A + w(^- 



9.9. 



where we define A = log (cosh + ^ sinh for simphcity of notation. For large values of TV 
as long as the ratio ^ and is finite, the behavior of the bound state energy is given by 

2 

(82) 

This problem again is an example of a certain kind of dimensional transmutation in non- 
relativistic quantum mechanics. By dimensional analysis, the hamiltonian of the system con- 
tains intrinsic energy scales and ^^^2 • However, after the renormalization, we obtain new 

parameters fif with energy dimensions. Hence, the number of parameters we expect for the 
energy at the beginning has changed after the renormalization. As it happens in the §^ case, 
the delta potentials on IP is an example of a generalized dimensional transmutation. 



Egr ~ 



log — log log 
log (cosh |r + I sinh I) 



5 Conclusion 

In this work, we studied a particle moving under the influence of A^ attractive Dirac delta 
interactions on some special Riemannian manifolds. We renormalized the problem and find a 
finite dimensional matrix <J>, called the characteristic matrix, by means of which a well defined 
expression for the resolvent can be written. All the information about the bound states can be 
obtained from the characteristic matrix. The renormalization can be done by means of the heat 
kernel and this is equivalent to the sharp cut-off method for the eigenvalues of the Laplacian, 
in the case of compact manifolds. We have studied the problem on particular compact and 
non-compact manifolds, S^, H^, and and we give explicit lower bound estimates on the 
bound state energies for each problem. Although we are concerned with particular manifolds, 
the basic idea for the renormalization can be applied also to general manifolds. 
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